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Abstract 

An important question in mathematical relativity theory is that of 
the nature of spacetime singularities. The equations of general relativ- 
ity, the Einstein equations, are essentially hyperbolic in nature and the 
study of spacetime singularities is naturally related to blow-up phenom- 
ena for nonlinear hyperbolic systems. These connections are explained 
and recent progress in applying the theory of hyperbolic equations in 
this field is presented. A direction which has turned out to be fruitful is 
that of constructing large families of solutions of the Einstein equations 
with singularities of a simple type by solving singular hyperbolic systems. 
Heuristic considerations indicate, however, that the generic case will be 
much more complicated and require different techniques. 

1 Introduction 

Anyone who has worked with nonlinear hyperbolic equations is familiar with 
the phenomenon that smooth solutions often cease to exist after a finite time. 
There are two common ways in which this happens. In the book of Alinhac|]lJ 
these have been called the 'ODE blow-up mechanism' and the 'geometric blow- 
up mechanism'. In the first case the solution itself is unbounded as the blow-up 
time is approached. In the second case the solution itself remains bounded 
while the first derivative is unbounded. Although the Einstein equations, which 
are the subject of the following, are arguably the most geometric of hyperbolic 
differential equations the blow-up which is typically encountered occurs by the 
ODE mechanism in the terminology of At the same time, the question of 
the sense in which the Einstein equations can be called hyperbolic is delicate. 
These points are related, as will now be explained 

The basic unknown in the Einstein equations is a pseudoriemannian metric 
g. The equations are diffeomorphism invariant, in the sense that if a metric g 
on a manifold M satisfies the Einstein equations and if <f> is a diffeomorphism 
of M onto itself, then the pull-back <p*g is also a solution. In order to make 
contact with PDE theory as it is usually formulated, it is necessary to choose a 
local coordinate system x a on M and treat the components g a p of the metric in 



this coordinate system as the basic unknowns. The Einstein equations imply a 
system of second order partial differential equations for these components. This 
system cannot have a well-posed Cauchy problem because of diffeomorphism 
invariance. For suppose that g a p is a solution of these coordinate equations 
with given initial data. The functions g a p are the components of a metric g. 
Now let 4> be a diffeomorphism which leaves a neighbourhood of the Cauchy 
surface invariant and let g = <jfg. Then g a p is a solution with the same initial 
data as g a j3. However it is in general not equal to g a p. Thus uniqueness in the 
Cauchy problem in the usual sense fails for this system. On the other hand, 
there is a different uniqueness statement which does hold. It says that if g 
and g are metrics with the same initial data on a suitable Cauchy surface then 
there exists a diffeomorphism <fi such that g — 4>*g which is the identity on the 
Cauchy surface. In the following the precise mathematical formulation of this 
statement will not be required, but the intuitive idea is important. It should 
be noted that this uniqueness property, often called geometric uniqueness, is 
exactly what is desired from the point of view of physics since metrics related 
by a diffeomorphism should be regarded as physically indistinguishable. 

It has now been indicated that diffeomorphism invariance presents a dif- 
ficulty when studying the Cauchy problem for the Einstein equations. How 
can this difficulty be overcome? In principle, one could attempt to construct a 
diffeomorphism-invariant version of PDE theory but this has apparently never 
been done. In practice, the procedure is to require some condition relating the 
coordinates to the metric. This breaks the diffeomorphism invariance and opens 
up the possibility that with this extra condition imposed the Einstein equations 
imply a system of hyperbolic equations called the reduced equations, and that 
conversely suitable solutions of the reduced equations give rise to solutions of 
the Einstein equations. Details on this procedure of hyperbolic reduction can 
be found in || . There are many different kinds of hyperbolic reduction and the 
task is to find reductions which are well-adapted to given problems. 

The Einstein equations are the basic equations of general relativity, which is 
the best existing theory of the gravitational field. Within this theory spacetime 
singularities, such as the big bang, correspond to singularities of g a p which 
cannot be removed by a diffeomorphism, even one which itself becomes singular 
where g a p does. 



A useful laboratory for studying questions concerning the global behaviour of 
solutions of the Einstein equations is the class of Gowdy solutions. These have 
high symmetry and after a suitable hyperbolic reduction give the following sys- 
tem of semilinear wave equations: 



2 The Gowdy equations 
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Here (t, x) £ (0, oo) x R and X and Z are real- valued functions. Because of the 
symmetry only one space dimension remains. It will be useful to compare this 
with the equation 

- d 2 u + d 2 x u = -e u (2) 

The equation (|J) and its higher dimensional analogues have been studied in a 
series of papers by Kichenassamy and LittmanjQ, [§[ |9|, |To| , [ll"| . Partial results of 
a similar type for ([!]) have been obtained in |l2[| and [|l5|. An obvious difference 
between the equations (|l|) and ^ is that while the singularities in solutions of 
(||) develop dynamically in the (t, x) plane there is a singularity in the coefficients 
of ([!]). Moreover it was proved by Moncrief[|l3j that smooth Cauchy data for (|j) 
at t = to > develop into a smooth solution on the whole of the interval (0, oo). 
In fact this difference is more apparent than real. When proving theorems 
about (||) it is convenient to introduce a function which is constant on the 
singular hypersurface and rewrite the equations in terms of that. In reality, a 
corresponding procedure has been carried out when introducing the coordinates 
used in the Gowdy equation. This is so natural in the context of the Einstein 
equations, where a priori all time coordinates are equally valid, that it might 
pass unnoticed. 

The nature of all these results is that they demonstrate the existence of 
large families of solutions of the equations whose singularities can be described 
precisely. This is achieved by using certain singular hyperbolic equations, whose 
singularity is of the type known as Fuchsian. The description of the singularity 
is given by an asymptotic expansion, which contains a number of free functions. 
These free functions may be called 'data on the singularity'. If these data are 
analytic then the expansion actually converges. The case of analytic data is 
easiest to handle. This was done for (|) in and § and for (Q) in The 
results for (||) were generalized to the case of smooth data (in fact data with 
sufficiently high order Sobolev regularity) in ||. The theory developed there 
does not suffice to cover (Q) with data which are smooth but not analytic. This 
difficulty was overcome in [fL5[ , where there is also a general discussion of the 
issue of extending results using Fuchsian techniques from the analytic to the 
smooth case and a survey of the literature where these techniques have been 
applied to problems in general relativity. 

Now the asymptotic expansions for solutions of the Gowdy equations will be 
presented for illustration. 

X(t,x) = X (x) +t 2k(x \^(x) +v(t,x)) (3) 
Z(t,x) = ft (a) log* + (f)(x) + u(t,x) (4) 

Here the functions u(t, x) and v(t, x) are remainder terms which are o{t) as 
t — > 0. Note the occurrence of powers of t with exponents depending on x in the 
expansion for X. It is these which cause additional technical difficulties in this 
case as compared to the case of (^]). Solutions of the above form are obtained 
for smooth data Xq, ft, <f> and ip which are arbitrary except for the restrictions 
that < k(x) < 1. The significance of the condition k > will not be discussed 
here. The condition ft < 1, known as the low velocity condition, will play a 
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significant role in the following discussion. It turns out that it can be removed 
at the expense of requiring that Xo should be constant. In this way a class of 
high velocity solutions can be obtained. They depend, however, on one less free 
function than the low velocity solutions. 

For the physical interpretation of the results on singularities in Gowdy space- 
times it is important to know that the singularity in the metric coefficients at 
t = cannot be removed by a diffeomorphism. One way of doing this is to 
consider scalar polynomial quantities in the curvature of the metric. These are 
real-valued functions which are defined in a coordinate invariant way. One of 
the simplest is the Kretschmann scalar K, which is the squared length (with 
respect to the given metric) of the Riemann curvature tensor. It is a function of 
X, Z and their first and second order partial derivatives whose explicit form is 
very complicated. Fortunately for a solution with an asymptotic expansion of 
the type given above it is not too difficult to calculate the leading term in the 
expansion of K about t = 0. The result, discussed in fl^| , is 

K ~ Q(fc)r 3 ( fc2+1 )/ 4 (5) 

for a certain quartic polynomial Q whose only positive root is at k = 1. Thus, 
provided k ^ 1 everywhere, K blows up uniformly as t — > 0. 

An important issue is that of the generality of the solutions constructed by 
the above techniques. They are general in the crude sense that they depend on 
the same number of free functions as data for the same equations on a regular 
Cauchy surface. In the case of (||) Kichenassamy proved in jl(| using the Nash- 
Moser theorem that a non-empty open set of data on a regular Cauchy surface 
develop into solutions of the type constructed by Fuchsian techniques. Up to 
now no corresponding results have been proved in problems arising in general 
relativity, such as the Gowdy problem. By a more direct approach, Chrusciel [|5| 
was able to show that a large open set of data for the Gowdy equations lead to 
solutions where the Kretschmann scalar blows up as t — > 0. Note that while the 
methods used in || were special to one space dimension, Fuchsian techniques 
do not depend essentially on the dimension. 

3 The Einstein-scalar field equations 

As mentioned above, the unknown of central importance in the Einstein equa- 
tions is a metric g. Physically it contains the information about the gravi- 
tational field and hence about the mutual gravitational attraction of material 
bodies. Thus it is natural that in any model involving the Einstein equations 
a description of the matter present is also necessary. In fact, in general rela- 
tivity the gravitational field has degrees of freedom which can propagate in the 
absence of matter. This leads to the phenomenon of gravitational waves. The 
equations in the absence of matter are known as the vacuum Einstein equations 
and have the simple geometric form Ric(g) = 0, where Ric(g) is the Ricci cur- 
vature of g. In the presence of matter these equations acquire a non-vanishing 
right hand side. An example which will be important in the following is that 
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of the Einstein-scalar field equations. In that case the Einstein equations take 
the form Ric(g) = V<^> ® V<f>. In general the matter fields must satisfy some 
equations of motion and in the case of the scalar field the equation of motion 
is that (f> should satisfy the linear wave equation defined by the metric g. In 
coordinate components this takes the form 

g a \d a d^ - ripd^) = o (6) 

where 

T l/3 = ^9 l5 (9 a g P s + d fj g aS - d 5 g a p) (7) 

are the Christoffel symbols of g. Here we use the summation convention and 
g a @ is the inverse matrix of g a p. 

It is known that in some cases the singularities of solutions of the Einstein 
equations display complicated oscillatory behaviour and it is suspected that this 
occurs under very general circumstances. The best known example of compli- 
cated behaviour is the Mixmaster solution. This is a solution of the vacuum 
Einstein equations which is symmetric under the group SU(2). The Einstein 
equations reduce to ordinary differential equations but the ODE solutions are 
very complicated. Until recently our knowledge of the Mixmaster solution was 
based on numerical calculations and heuristic arguments. Now a number of 
features of this picture have been proved rigorously by Ringstr6m[[l6]. The gen- 
erality of Mixmaster behaviour is still out of reach of present analytical results. 
There have, however, been interesting advances in numerical approaches. (See 

General predictions about the nature of singularities in solutions of the 
Einstein equations come from heuristic work of Belinskii, Khalatnikov and 
Lifshitzj3j. According to their picture the singularities in solutions of the vac- 
uum Einstein equations should look like a different Mixmaster solution at each 
spatial point near the singularity. In particular the evolution at different spa- 
tial points decouples and the solution of the PDE system can be approximated 
by solutions of an ODE system. (Here we are reminded of the ODE blow-up 
mechanism in the terminology of Alinhac.) The same picture suggests that for 
the Einstein-scalar field system the decoupling should still take place but the 
ODE solutions should now be monotone rather than oscillatory near the sin- 
gularity. This gives rise to a hope that the Einstein-scalar field system should 
be easier to handle analytically than the Einstein vacuum equations. Following 
this lead, it was shown in Q that a Fuchsian analysis could be carried out for 
the Einstein-scalar field system without any need for symmetry assumptions. 
Solutions are obtained which depend on the same number of free functions as 
there are in the Cauchy data on a regular Cauchy surface. In this crude sense 
the solutions are general. Only the analytic case was treated up to now. 

The asymptotic expansions obtained will now be presented. This is only 
done for a certain diagonalizable case. This does not restrict the number of free 
functions although it is a restriction on the solutions and a large part of the 
work of [Ej was devoted to showing that the restriction is unnecessary. On the 
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other hand, the diagonalizable case is best suited to explain the results. For 
physical reasons we are interested in solutions on a four-dimensional manifold. 
The time coordinate is denoted by x° and the spatial coordinates by x a with 
a = 1, 2, 3. The Gaussian coordinate conditions 

5oo = -1 (8) 
9oa = (9) 

are imposed. The asymptotic expansions for the remaining unknowns are as 
follows: 

gab = t 2 PH a l b + t 2 P 2 m a m b + t 2 P3n a n b + o(t q ° b ) (10) 
4> = A\ogt + B + o{\) (11) 

Here pi, P2, P3, l a , m a, n a , A and B depend on the spatial coordinates x a and 
Pi +f>2 +P3 = 1- The q ab are suitably chosen real numbers, which will not be 
discussed in detail here. The important thing is that they represent higher order 
terms. The functions occurring in the expansion must satisfy certain algebraic 
and differential equations which will also not be discussed, except to say that 
there is a good understanding of the solution set of these equations. There is an 
equivalent of the condition k < 1 in the Gowdy case, namely that all p a should 
be positive. This implies that all p a are smaller than one. 

The conclusion of this work is that it provides strong support for the ideas 
of Bclinskii, Khalatnikov and Lifshitz in the case that a scalar field is present. 
The conceptual, geometric and calculational aspects of the proof are a lot more 
complicated than in the Gowdy case, but the analytic input is almost identical. 



4 The nonlinear scalar field 

From the foregoing discussion of the great difference in the nature of singularities 
in solutions of the Einstein equations caused by the presence of a scalar field 
the impression might arise that the whole question depends in a delicate way on 
the exact nature of the matter fields present. In fact this impression is probably 
false: there should probably only be a small number of universality classes which 
cover most types of matter. This will be illustrated by a class of matter models 
which will be referred to collectively as the nonlinear scalar field. It turns out 
that a wide variety of choices of the parameters which go into the definition all 
lead to the same asymptotic behaviour, which is that of the (linear) scalar field 
discussed in the previous section. This will also allow some of the aspects of the 
Fuchsian theory to be shown in action. 

The main interest for physicists of the equations obtained by coupling the 
Einstein equations to a nonlinear scalar field has been the the phenomenon 
known as inflation. The mathematics of relevance to inflation concerns the 
behaviour of the solutions in the direction away from the singularity and so has 
no direct connection with the discussion below. 
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Consider first the following equation 

-d 2 t 4>-r 1 d t (t> + Ac/) = m 2 (f> + V'{(t)) (12) 

Here (pis a, real-valued function, A is the Laplacian on R n (the case of n = 3 is of 
most interest for the following) and to is a constant (mass) which we may as well 
take to be non-negative. The function V (potential) is smooth and V vanishes 
at least quadratically at the origin. The separation of the right hand side into 
the sum of a mass term and a potential term has no essential significance in 
this problem. It is convenient for making contact with the notation used in the 
physics literature. The aim is to determine whether this equation has solutions 
which near t = have the leading order asymptotics A(x) logi + B(x) already 
encountered in the previous section. To this end, introduce ip = t~ e (<p — Alogt — 
B), where e is a positive constant, whose purpose will be seen later on. When 
reexpressed in terms of ip the equation jl^ ) becomes: 

- t 2 d?tp - (2e + l)td t tp -e 2 tp + t 2 Aifj = -t 2 ~ e log tAA - t 2 ~ e AB 

+t 2 - e m 2 {A logt + B + t c ip) + W(t, x, i/0 (13) 

where W(t,x,ip) = t 2 ~ e V'(A\ogt + B + t e ip). The equation has been multiplied 
by t 2 so as to make the combination td/dt appear. Some choices of V to be found 
in the physics literature are V\(<p) = Vq, a constant, V2{<p) — 4> a , Vz{<j>) = 
for positive constants a and j3. The value of the constant Vq has no effect on 
the equation (|l^) but does have an effect when the coupling to the Einstein 
equations is considered. The potential V\ combined with a non-zero value of to 
occurs in a model known as chaotic inflation. The potential V3 occurs in what 
is called power-law inflation. (It is not the potential itself which has power- 
law behaviour, but the long-time behaviour of the solution.) The functions W 
corresponding to V\, V% and V3 are W\{t, x, ip) = 0, Wi(t, x, ip) = at 2 ~ e (A\ogt+ 
B + Fip)™- 1 and W 3 (t,x,ip) = [3e tm t 2 -' i+l3A exp(/3t e V)- The functions W u W 2 
and W3 are smooth for t > 0. Suppose that the function A is nowhere vanishing. 
Then W% and Wi extend continuously to t — by defining them to be zero there. 
This is also true of their derivatives with respect to the arguments x and ip. This 
will be abbreviated by saying that the functions W\ and W2 are regular. As for 
W3, can be made regular in this sense by choosing e sufficiently small provided 
A > —2(3~ 1 . This is a restriction on the data on the singularity. In particular, 
if (3 is negative then it bounds A awav from zero. If V had been chosen to grow 
even faster at infinity, e.g. V(<fi) — then the corresponding W could not have 
been made regular by imposing inequalities on A. This type of potential with 
faster then exponential growth does not seem to occur in the physics literature 
and the Fuchsian theory does not seem to apply to it. It will not be considered 
further here. 

The Fuchsian theory we will use concerns first order systems and thus the 
first step in the analysis is to reduce ([l3]) to first order in a suitable way. Let 
v = tdtip + &p and w = td x ip. Then the the following system is obtained: 

tdti> + eip — v = 
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td t w + (e - l)w = td x v (14) 
tdtv + ev = td x w + AAt 2 ~ e logt + t 2 ~ £ AB 

-t 2 - e m 2 {A\ogt + B + i/j) + W(t,x,ip) 

This system is of the following general form for Fuchsian equations: 

tdtU + N(x)u = t r if{t,x lU ,d x u) (15) 

for some rj > 0. It is also symmetric hyperbolic, a luxury which is usually not 
easy to obtain in problems in this area. On the other hand, it does have a 



disadvantage. The theorem proved in 12 required the assumption that 



be bounded for t near zero. To obtain this the eigenvalues of N(x) should have 
non-negative real parts. Since e has to be chosen small in general, depending 
on the data, the eigenvalue e — 1 will be negative and this condition is not 
satisfied for the equation (|lj). An alternative form is obtained by substituting 
the definition of v back into (|TJ) in one place. The result is: 

td t ifj + etfj — v = 

td t w + ew — td x v + td x if) (16) 
td t v + ev = td x w + AAt 2 ~ e logt + t 2 ~ e AB 

—t 2 ~ e m 2 (Alogt + B + ip) + W(t, x, ip) 

The equation ( [l6| ) is not symmetric hyperbolic. It does, however satisfy the 
conditions of the existence theorem for the case of analytic data proved in Q . 
In order to ensure this it is important that e is positive. If e were chosen to be 
zero then all the eigenvalues of the matrix N(x) would vanish but the matrix 
would not be diagonalizable. This would violate one of the conditions assumed 
in [ ]l2[ , namely that zero eigenvalues, if they occur at all, only occur through 
N(x) being a direct sum of a matrix whose eigenvalues have positive real parts 
with a zero matrix. 

The result of the previous discussion is that given analytic functions A and B 
with A nowhere vanishing then if the restrictions on V and A already obtained 
are satisfied there exists a unique solution, analytic in x and continuous in t 
of ([l2|) with the desired asymptotics. This shows the existence of a family of 
solutions depending on two free analytic functions, which is the same as the 
number of functions which can be given as data on a regular Cauchy surface. 

In []l5) a method for going beyond the analytic case in Fuchsian problems 
was applied to the Gowdy equations. Since the matrix N(x) occurring in ( [l6| ) 
does not depend on x it should be possible to prove an existence theorem for 
that equation with smooth data given on the singularity by the methods of Q . 
Here we would like to note that the method of |l5| generalizes easily to the case 
of smooth data for (|l^). The main idea is to use both the systems ( fLl| ) and 
(^6|). If the smooth data is approximated by a sequence of analytic data then a 
corresponding sequence of analytic solutions can be produced using (|l|). These 
define a sequence of analytic solutions of (jl(]). Then the fact that the latter 
equation is symmetric hyperbolic allows the use of energy estimates to prove 
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that the sequence of analytic solutions converges to the desired smooth solution. 
Details of the procedure can be found in jl5| . 

5 Coupling the nonlinear scalar field to the Ein- 
stein equations 

In this section it will be shown that the Einstein equations coupled to the 
nonlinear scalar field can be handled by the same method as that used for the 
case with m = and V = in |J. The arguments will only be sketched since 
most of the steps are identical. The potential will be assumed to be one of the 
functions Vi, V2 or V3 considered in the last section. In addition it is assumed 
that A 7^ and that in the case of V3 the inequality A > — holds. The 
Einstein equations take the form 

Ric( 5 ) = V<£ ® V<t> + [^rn 2 (f> 2 + V(<p)]g (17) 

Evidently this equation does change if V is changed by an additive constant. 
The equation for the scalar field is 

g a0 (d a d f ^ - Vl p d l( j>) = m 2 + V'(ct>) (18) 

Applying the Gaussian coordinate conditions 500 = — 1 and go a = gives 

g al3 d a d cf> = -dU + g ab d a d b <p (19) 

and the close relation of equation ( |l8| ) to equation ([l2]) begins to become ap- 
parent. We are looking for solutions to the Einstein equations coupled to a 
nonlinear scalar field where the metric has the asymptotic form found in 
for the case of a linear scalar field. This form implies that the quantities t 2 g ab 
all vanish like some positive power of t as t — > and thus do not disturb the 
Fuchsian form of the wave equation. The same holds for the quantities tT a and 
the quantity tT° is equal to one up to higher order terms. These facts cannot be 
easily seen from the equations given in this paper. To check them it is necessary 
to use more detailed results from Q . The conclusion is that up to higher order 
terms the equation (|l^ ) agrees with the model equation JI^). 

What remains to be done in order to check that the Einstein equations 
coupled to a nonlinear scalar field have solutions which behave asymptotically 
like those for a linear scalar field is to see that certain source terms in the 
Einstein equations are equal to those in the case with a linear massless scalar 
field (which has already been solved in up to terms of higher order in t. 
The essential thing is to estimate the quantity d 2 </> + ^m 2 (f> 2 + V((f>). More 
specifically it is necessary to estimate the deviation of this quantity from the 
corresponding quantity with -0 = and show that it is O(t~ 2+ao ) where ag 
is a positive parameter introduced in In the course of the proofs in that 
paper it is permissible to reduce the size of the number a$ if required. It 
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is easy to check that, under the assumptions already made on V and A, the 
parameter cxq can be chosen so as to ensure this estimate. This is the same kind 
of procedure as choosing e small in the case of the nonlinear scalar held. The 
information obtained in this way must be used in two places. It must be used 
in the Einstein constraints in estimating the quantity C of 0. It must also be 
used in evolution equations for in 0, which constitute the essential part 
of the Einstein evolution equations. The additional contribution which arises 
when passing from the massless linear scalar held to the nonlinear scalar held 
does not change the equations for a / b. For a — b the coefficient ao appears. It 
can be concluded from all this that solutions of the Einstein equations coupled 
to a nonlinear scalar held which have the same asymptotics as in the case of 
a massless linear scalar held and which are of the same degree of generality 
can be constructed by Fuchsian techniques. This has been shown under the 
assumption that the potential is of one of the forms Vi , Vi or V3 and that in the 
third case the initial datum A satisfies an additional restriction. In this proof 
the velocity-dominated system, which is the system of equations satisfied by the 
leading order approximation to the solution, is the same as in the case of the 
massless linear scalar held. 



6 Formation of localized structure 

It has already been indicated that the use of Fuchsian theory to construct space- 
times with singularities of a well-understood type is not likely to be sufficient 
to handle general singularities of solutions of the Einstein equations, due to the 
appearance of oscillatory behaviour in time as the singularity is approached. 
There is, however, another effect which complicates the situation in an essential 
way. Note that the solutions constructed by Fuchsian techniques are such that 
the solution blows up like log t as t = is approached and the spatial derivatives 
of the solutions of all orders are also O(logi). Taking spatial derivatives does 
not increase the rate of blow-up. It will now be shown that there are solutions 
for which this property does not hold. Consider the transformation: 

X = X/(X 2 + e 2Z ) (20) 
Z = log[e z /(X 2 + e 2Z )} (21) 

This transformation is equal to its own inverse and transforms solutions of the 
Gowdy equations into solutions. The idea is now to take a solution {X, Z) with 
the low velocity asymptotic behaviour and transform it to a new solution [X, Z) . 

The case will be considered that Xq(0) — 0, Xq(x) ^ for x 7^ and 
d x X$(Q) ^ 0. For the transformed solution X 2 + e 2Z is of the form 

(X ) 2 + (e 2 * + 2X j>)1?~ k + o{t 2k ) (22) 

From this it is not hard to see that X is bounded on an interval about x = 
uniformly in t. Similarly Z = O(logi). Thus the transformed solution (X,Z) 
does not blow up any faster than the original solution (X, Z). However spatial 
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derivatives blow up faster in general. To see this it suffices to evaluate Z(t, 0) 
and Z(M fe/2 )- The second of these is bounded for values of x close to x = 
as t — > while the first behaves like — klogt. It follows from the mean value 
theorem that the maximum value of d x Z in a small region about x — must 
blow up at least as fast as t~ k l 2 . Of course precise rates of blow-up could be 
calculated for this and higher order derivatives if desired. It is seen that these 
solutions of the Gowdy equations, unlike those constructed by the Fuchsian 
method, have spatial derivatives which blow up faster than the function itself. 
There is formation of localized structure. 

In the explicit example given above the localized structures formed do not 
have an invariant geometrical meaning and are due to a bad choice of variables. 
On the other hand, there is convincing numerical evidence due to Berger, Mon- 
crief and collaborators for a second kind of localized structure. In the latter 
case there is no reason to doubt that there is a real geometrical effect involved. 
For a survey of this work and some pictures of the structures, see [|J. There is 
also a lot of interesting work on the subject using heuristic methods. A rigorous 
treatment is still lacking. 

Although this is still in the domain of speculation, there are arguments that 
in more general situations things will get yet more complicated. The basic 
idea will now be sketched. The oscillatory behaviour in the Mixmaster model is 
composed of a series of bounces where expansion turns into contraction and vice 
versa. In the inhomogeneous case bounces of this kind are supposed to happen 
independently at different spatial points. However this does not always happen 
coherently. There are spatial regions of coherent behaviour but these split into 
smaller regions as the singularity is approached. On the boundaries between 
regions spatial derivatives become large. The localized features observed in the 
Gowdy solutions are simple examples of boundaries of this kind. The heuristic 
models available suggest that in the Gowdy case this kind of formation of struc- 
ture will only happen a finite number of times, at least for generic solutions. In 
the general case, on the other hand, it is expected to happen infinitely many 
times before the singularity is reached, creating structures on arbitrarily small 
spatial scales. This behaviour has been called 'spacetime turbulence'. Informa- 
tion on this can be found in j|3]. It has also been suggested that this process 
which is predicted by the picture of Bclinskii, Khalatnikov and Lifshitz actually 
shows that the picture must eventually break down in the generic case before 
the singularity is reached. 

It seems clear that it will take a long time before the ideas mentioned in the 
last paragraph can be brought into the fold of well-understood mathematical 
phenomena. Nevertheless, the recent progress described in this paper is a ground 
for optimism that there will be significant mathematical developments in the 
area of singularities of solutions of the Einstein equations in the next few years. 

References 



li 



Alinhac, S. 1995 Blowup for nonlinear hyperbolic equations. Birkhauser, 
Boston. 

Andersson, L., Rendall, A. D. 2000 Quiescent cosmological singularities. 



Preprint gr-qc/0001047 



Belinskii, V. A., Khalatnikov, I. M. and Lifshitz, E. M. 1982 A general 
solution of the Einstein equations with a time singularity. Adv. Phys. 31, 
639-667. 

Berger, B. K. 1997 Numerical investigation of singularities. In M. Fran- 
caviglia et. al. (eds.) Proceedings of the 14th International Conference on 
General Relativity and Gravitation. World Scientific, Singapore. 

Chrusciel, P. 1991 On uniqueness in the large of solutions of the Einstein 
equations ('Strong Cosmic Censorship') Proc. CMA 27, ANU, Canberra. 

Friedrich, H., Rendall, A. D. 2000 The Cauchy problem for the Einstein 
equations. In B. G. Schmidt (ed) Einstein's Field Equations and Their 
Physical Implications. Lecture Notes in Physics 540. Springer, Berlin. 

Kichenassamy, S. and Littman, W. 1993 Blow-up surfaces for nonlinear 
wave equations I. Comm. PDE 18, 431-452. 

Kichenassamy, S. and Littman, W. 1993 Blow-up surfaces for nonlinear 
wave equations II. Comm. PDE 18, 1869-1899. 

Kichenassamy, S. 1996 Fuchsian equations in Sobolev spaces and blow-up. 
J. Diff. Eq. 125, 299-327. 

Kichenassamy, S. 1996 The blow-up problem for exponential nonlinearities. 
Comm. PDE 21, 125-162. 

Kichenassamy, S. 1996 Nonlinear Wave Equations. Marcel Dekker, New 
York. 

Kichenassamy, S., Rendall, A. D. 1998 Analytic description of singularities 
in Gowdy spacetimes. Class. Quantum Grav. 15, 1339-1355. 

Moncrief, V. 1980 Global properties of Gowdy spacetimes with T 3 x R 
topology. Ann. Phys. (NY) 132, 87-107. 

Montani, G. 1995 On the behaviour of the universe near the cosmological 
singularity. Class. Quantum Grav. 12, 2505-2517. 

Rendall, A. D. 2000 Fuchsian analysis of singularities in Gowdy spacetimes 



beyond analyticity. Preprint gr-qc/0004044 



Ringstrom, H. 2000 Curvature blow-up in Bianchi VIII and IX vacuum 
spacetimes. Class. Quantum Grav. 17, 713-731. 

Weaver, M., Isenberg, J. and Berger, B. K. 1998 Mixmaster behavior in 
inhomogeneous cosmological spacetimes. Phys. Rev. Lett. 80, 2984-2987. 



12 



